home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dpptri.z / dpptri
Encoding:
Text File  |  2002-10-03  |  3.4 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPPPPPTTTTRRRRIIII((((3333SSSS))))                                                          DDDDPPPPPPPPTTTTRRRRIIII((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPPTRI - compute the inverse of a real symmetric positive definite matrix
  10.      A using the Cholesky factorization A = U**T*U or A = L*L**T computed by
  11.      DPPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DPPTRI( UPLO, N, AP, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, N
  19.  
  20.          DOUBLE         PRECISION AP( * )
  21.  
  22. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  23.      These routines are part of the SCSL Scientific Library and can be loaded
  24.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  25.      directs the linker to use the multi-processor version of the library.
  26.  
  27.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  28.      4 bytes (32 bits). Another version of SCSL is available in which integers
  29.      are 8 bytes (64 bits).  This version allows the user access to larger
  30.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  31.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  32.      only one of the two versions; 4-byte integer and 8-byte integer library
  33.      calls cannot be mixed.
  34.  
  35. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  36.      DPPTRI computes the inverse of a real symmetric positive definite matrix
  37.      A using the Cholesky factorization A = U**T*U or A = L*L**T computed by
  38.      DPPTRF.
  39.  
  40. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  41.      UPLO    (input) CHARACTER*1
  42.              = 'U':  Upper triangular factor is stored in AP;
  43.              = 'L':  Lower triangular factor is stored in AP.
  44.  
  45.      N       (input) INTEGER
  46.              The order of the matrix A.  N >= 0.
  47.  
  48.      AP      (input/output) DOUBLE PRECISION array, dimension (N*(N+1)/2)
  49.              On entry, the triangular factor U or L from the Cholesky
  50.              factorization A = U**T*U or A = L*L**T, packed columnwise as a
  51.              linear array.  The j-th column of U or L is stored in the array
  52.              AP as follows:  if UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for
  53.              1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for
  54.              j<=i<=n.
  55.  
  56.              On exit, the upper or lower triangle of the (symmetric) inverse
  57.              of A, overwriting the input factor U or L.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDPPPPPPPPTTTTRRRRIIII((((3333SSSS))))                                                          DDDDPPPPPPPPTTTTRRRRIIII((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      INFO    (output) INTEGER
  75.              = 0:  successful exit
  76.              < 0:  if INFO = -i, the i-th argument had an illegal value
  77.              > 0:  if INFO = i, the (i,i) element of the factor U or L is
  78.              zero, and the inverse could not be computed.
  79.  
  80. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  81.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  82.  
  83.      This man page is available only online.
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.